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Abstract: Using a technique, analogous to coherent population trapping 
in an atomic system, we propose schemes to create transverse light propa- 
gation violating left-right symmetry in a photonic circuit consisting of three 
coupled waveguides. The frequency windows for the symmetry breaking of 
the left-right energy flow span over 200 nm. Our proposed system only uses 
linear passive optical materials and is easy to integrate on a chip. 
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1. Introduction 

The integration of nonreciprocal photonic devices on Si or CMOS platforms has been chal- 
lenging in the past decades. The known optical nonreciprocity can be divided into two classes: 
forward-backward nonrecipricity (FBNR) and left-right nonrecipricity (LRNR). The first class 
of nonreciprocal component can be realized by various approaches using the magneto-optic ef- 
fects |[l]|2l|3]|4l|5]|6|, nonhnearity Q, a modulation media EHHITOl, optomechanics lHHIHl, or 
magnetized plasmonic metal lfT3l[T4l . It can be used for optical isolators. To date, the LRNR, 
indicating the nonreciprocal light flow between left and right ports of photonic circuits, has 
been discussed in an array of coupled waveguides by only a few research groups ifTS] [T6] [TtI . 
If the complex optical potential causes the parity-time (PT) symmetry breaking, two coupled 
waveguides can show a LRNR light transfer in the transverse direction llT5l[T6l . The realization 
relies on the precise control of the active medium. Although unclear for optical isolators, the 
photonic circuit with the left-right nonrecipricity may merge the incident light beams. 

The need for integration of optical nonreciprocal elements on a Si/CMOS chip platform is a 
long-standing problem. The realization of the nonreciprocal light propagation in a completely 
linear optical medium can strongly impact on both fundamental physics, and also vast applica- 
tions for integrated optics because of the compatibility with the Si material and CMOS chips. 
Again based on the PT symmetry breaking induced by a periodic modulation of complex op- 
tical potential, Feng et al. stated that they, for the first time, observed the nonreciprocal light 
propagation in a linear passive optical material llTSl . Unfortunately, they admitted their mistake 
|fT9l after Fan et al. commented on their work Il20l . Another experimental realization of on-chip 



optical diodes using all-dielectric, passive, and linear silicon photonic crystal structures is re- 
ported by Wang et al. 112111221 . They confidentially explained why it is possible to make optical 
diodes using a spatial symmetry breaking geometry in a passive and linear optical medium ll22l . 

In the LRNR, the input light is always localized in one waveguide ifTSl [161 [TTl . This is 
similar to the coherent population trapping in a three-level A-type atom Il23l . While the bend- 
ing waveguide array can simulate well the quantum dynamics of atoms. The classical optical 
analogs of coherent population transfer ll24l l25l and population trapping in the continuum 1261 
and Rabi oscillation 1271 as well has been proved by Longhi's group. Although the trapping of 
equal light in two waveguides has been discussed l24ll . the results did not show a valid LRNR 
of light flow. A theory work indicated that the critical large nonlinearity is necessary to induce 
nonrecipricity ifTTl in two evanescently coupled waveguides. However the LRNR in the trans- 
verse light flow has been recently observed ifTS] [T6l . Moreover, the three/many -body systems 
behavior essentially different from the simple two-body system studied in [17]. At least, the 
optical trapping in coupled three waveguides analogous to the atomic CPT can not be achieved 
in an optical system composing of two waveguides. It is interesting if one can realize the non- 
reciprocal wave propagation in evanescently coupled linear and passive waveguides. We expect 
to achieve the LRNR in a waveguide array. This is the motivation of our work. 

Here we propose simple methods to generate the second class of optical nonreciprocity in 
an array of three coupled waveguides only making from linear, passive optical materials. We 
focus on the nonreciprocal transverse energy flow between left and right optical waveguides. 
Thanks to a small dispersion of a linear waveguide, our system can behave in a nonreciprocal 
manner in an ultrabroad band. The results by numerical simulation of beam propagating method 
(BPM) and solving the coupled mode equation (CME) demonstrate the breaking of symmetry 
of transverse energy flow. 

2. Setup and model 

Our system, shown in Fig. [1] is composed of three coupled waveguides embedded in a wafer 
of width W and length L. The middle waveguide D3 couples to the waveguides Di and D2. We 
assume that the coupling between the waveguides Di and D2 is negligible. We also assume that 
the two side waveguides are lossless but some loss can be included in the middle one. In our 
photonic system, eigenmodes in the individual waveguide exchange energy via their evanescent 
fields when two waveguides are close. The couplings are denoted as kij, and K23, and decrease 
as the distance t/i3,23 between two waveguides increases. The coupling between Di and D2 is 
assumed vanishing because these two waveguides are far enough from each other The field in 
D3 decays exponentially with a constant y that can be controlled lITSl [TSl [161 . 

The evolution of power of field in photonic circuits can be studied either by numerical solving 
the Helmholtz equation or with the derived coupled mode theory. However the later presents a 
clearer physic understanding. We first present the approach based on Helmholz equation. 

We consider our system to be two-dimensional (2D). This is reasonable if the size of waveg- 
uide in the y direction is much larger than that in the x direction. The propagation of the field E 
in the photonic circuits in a 2D space can be described by the Helmholtz equation, which takes 
the form 

d^E d-E , 

^ + ^+^o'^(-'^'^)£-0, (1) 

where z is the propagating direction, x is the transverse direction, and A:o = ^ is the wave vector 
of field with wavelength A in the free space. The dielectric constant e(x,z) plays the role of the 
optical potential. We can resolve the field E into its slowly varying amplitude 41 and a fast 
oscillating factor, E - ij/e'^^P' , where the propagation constant p - kQiigff. tieff is the effective 
index of waveguides. The sign before /3 indicates the propagating direction; minus (plus) for the 
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Fig. 1. Light trapping in a photonic circuit consisting of three waveguides embedded in a W 
wide, L long substrate. The coupling /ci3 « ko^. rfi3,23 is the distance between two waveg- 
uides. The dielectric constants of substrate is e^, while (; 6 j 1, 2, 3)) is the constant profile 
of corresponding individual waveguides D,- without others. The widths of waveguides are 
tj, respectively. 



propagation along the positive (negative) z direction. In the paraxial approximation |2 i/6^| » 

I d^il/ 1 

L the amplitude of the field in the photonic circuit evolves according to the equation 

We numerically solve Eq. (|2|i with BPM to simulate the propagation of field with a spatial 
resolution & = 1 ixm and 6x = 0.1 i_im. A finer spatial grid gives the same results. Throughout 
simulation, we use a Gaussian input profile i//{x,X(),z - 0,L) - exp(-(jc- jcq)^ /2wp), where 
Wp - 1 \im is the half width of waist of input field, xq = Xq' (/ 6 (1,2)) is the center position 
of input ports Pli.l2,r\- The profiles ip{x,x^^\Q) and \lr(x,x^^\L) corresponds to the input field 
launching into the waveguide Di from the port Pn at z = and Pr\ slz — L, respectively. While 
the profile ij/{x,x^ ,0) means an input to the port P12 at z = 0. The intensity of the field at the 
peak is unity. This profile is very close to the fundamental eigenmode of Di or D2. 

Our simulation focuses on a hght with wavelength Aq - 1.55 iJ.m , which is of interest 
in optical communications. The photonic circuit can be integrated in a wafer with a sub- 
strate dielectric constant Sj - 10.56eo, where eo is the permittivity of free space. We con- 
sider weakly guiding waveguides with ecme = 10.76eo in core. The imaginary part of the 
dielectric constant in the middle waveguide D3 is V>{e^] - -O.Oleo- This induces a loss of 
y - 1 1.6 mm~^ according to our numerical simulation. The propagation constant is calculated 
by solving the eigenvalue equation 1281 of TE mode. It is evaluated to be |y6| — 1.32593 x 10^ 
rad/m~ k^i yjEcre + y/sT)/^ yfsQ for all simulation. This is a good approximation for the weakly 
guiding waveguides. 

In our numerical method, the Neumann boundary condition (NBC) is used to greatly suppress 
the reflection field from the transverse boundary. A small reflection, which can be a practical 
noise from the boundary of device in experiments, is responsible for the background noise of 
our numerical results. 

Before discussing the results, we present the equivalent, but physically transparent, coupled 
mode equation approach to explain our system. 

The light field E in photonic circuits can be expressed as a supermode of eigenmodes £, of 



individual waveguide D,, i.e., 



Yj Z ^'-A/e^'^'S (3) 

i=l,23 i'=l,2,3 

where the amplitude of eigenmode in the rth waveguide is denoted by A,- (/ e {1,2,3)). Here tf/i is 
the slowly varying envelope of eigenmode £, and is the corresponding propagation constant, 
which can be controlled by designing the dielectric constant e,- of individual waveguide D, and 
its width f,. This parameter also depends on the wavelength of light. According to the coupled 
mode theory 1281 l29l l30l derived from the Helmholtz equation, the dynamics of the modal 
amplitudes is described by 

jAisAi +>i3(z)A3 (4a) 

jA23Ai+jK23{z)Ai (4b) 

>i3(z)Ai + jK23iz)A2 - yAi , (4c) 

where An = /3i -ySs and A23 = P2 -Pi are the phase mismatch between waveguides Di (D2) 
and D3 . The power in the /th waveguide is evaluated by Pi = A,A* . The total power in the system 
is P = YjiPi- The phase mismatching and coupling in the coupled mode equation Eq. (|4]i can be 
derived from Helmholtz equation Eq. These parameters are dependent on ko and the optical 
potential s{x,z)- The coupling is given by 1281 l29l l30l 

where s denotes the cross section of space and 6s„ - s{x,z)-s„ at the cut position z- For 
simplicity's sake, we have assumed weakly guiding waveguides and the relation k„,„ = k„,„. We 
also neglect the second-order spatial derivatives of the amplitude A, and the small self phase 
shifts due to the perturbation of neighbor waveguides. A full study of the relation of parameter 
to the Helmholtz equation has been presented by Hardy et al. l29l[30l . Note that the coupled 
mode theory presents a general model. In contrast, the numerical results depend on the structure 
of system, and only provides one of many implementations. The different structures can lead 
to the same set of parameters in the coupled mode equation. On the other hand, the numerical 
simulation presents a full picture of light in photonic circuits. 

Next we turn to our idea about how to create optical nonreciprocity in the transversal energy 
flow in three coupled waveguides by giving a connection of our classic photonic circuit to a 
quantum system. 

Due to the equivalence between the Helmholtz equation in photonic circuits and the 
Schrodinger equation in quantum mechanics, the behavior of light propagating in a photonic 
circuit is similar to the dynamics of the internal atomic states of a quantum system l24l l25l . 
For example, the normalized light power trapped in optical waveguides plays the role of atomic 
population. Optical nonreciprocities in an array of coupled waveguides can then be considered 
as the trapping of input light on demand. As is well known, in a A-type three-level atomic sys- 
tems, we can adiabatically create a target state independent of the initial state of system via the 
so-called coherent population trapping (CPT) ll23l . Our system is analogous to such a A-type 
three-level system. Just as in CPT in the atomic system, we expect to trap most light energy in a 
selected optical waveguide by suitably controlhng the coupling between the waveguides. This 
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is the basis of the optical nonreciprocity studied in this paper. Note that the LRNR we propose 
here is substantially different from the FBNR used for optical isolator on the basis of the break- 
ing of Lorentz reciprocity theorem. In the former case, both two sources input into ports in the 
left hand side but their responses coming out from the right hand side are compared, while the 
source and the response must exchange in the later Thus the implementation of the LRNR in a 
linear, passive medium is not forbidden by the Lorentz reciprocity theorem, which requires the 
exchange of the place of source and response. 

The structure of photonic circuit to create the nonreciproc al transverse energy flow of light 
is shown in Fig.|2](a). Through our system, we assume no loss in waveguide Di and D2. The 
dielectric constant is Ss = 10.56eo in the substrate, while it is Score = 10.76eo in core of wi 
and W2. In waveguide wj, the dielectric constant in core is sj = (10.76 - 0.01/)eo. The straight 
waveguide W3 is 4 |j.m wide around its center W3(0) - 26.7 |j.m. The waveguides wi and W2 
with width 2 |i.m are curves along their varying central positions, which are defined as 



wi(z) [|xm] 



31; z < 1 . 1 mm 

31+3(l+sin(2;r(z-1350)/1000))/2.0; z> 1.1 mm&z< 1.6 mm 
34; z > 1.6 mm 



and 



W2iz) [l^m] 



22; 



z < 0.1 mm 



22 + 1.4(1 +sin(27r(z-600)/2000))/2.0; z>0.1 mm&z< 1.1 mm 



23.4; 



z > 1 . 1 mm 



To fit the numerical results, the coupling ki^ and K2j, are assumed to vary corresponding to the 
central positions wi (z), W2(z) and W3(z) of waveguides. The other parameters for coupled mode 
equation are given by: 



y -10 mm , 
Ai3 =A23 = 20 mm" 



(6) 
(7) 



To simulate the varying gaps between waveguides, we assume two gradiantly changing cou- 
pling strength /C13 and a'23 as 



a:i3(z) [mm ] = 



and 



Kiiiz) [mm ] = 



25; z < 1.1 mm 

l+24(l-sin(27r(z-1350)/1000)/2.0; z> 1.1 mm&z< 1.6 mm , 
1 ; z > 1 .6 mm 



2; z < 0.1 mm 

2 + 88(1 +sin(27r(z-600)/2000))/2.0; z>0.1 mm&z< 1.1 mm 

90; z > 1 . 1 mm 



The corresponding coupling strength in the propagating direction for Eq.|4]are shown in Fig.|2] 
(b). Two waveguides in the same chip always couples to each other even if the coupling strength 
is very small. To consider this coupling, we assume small values as the distance between two 
waveguides are large. The intensity in waveguides given by the coupled mode equations change 
very slightly if we neglect this small coupling. We note that, in the absence of loss in the 




Fig. 2. (a) The waveguide structure for left-right nonrecipricity. (b) The coupling as a func- 
tion of propagating distance z. Blue solid lines (i) for /C13, red dashed lines (ii) for K23. Here 
Ao = 1.55 \im. 



waveguide W3, the system is reciprocal (not shown here). However if we include loss in the 
middle waveguide, we create left-right nonreciprocities. 

Let us assume that a light with unit amplitude is incident on the port Pn or f ^2 at z = 0. 
If the photonic circuit is reciprocal, the transmission from port Pri or Pr2 exchanges as well 
if the incident exchange. However, in the case of left-right nonrecipricity, the light launched 
into port Pn and P12 always effectively transfer to the waveguide Di and comes out from port 
Pri. If we use constant couplings ki^ and K23, the LRNR is obtained but the transmissions are 
small. The energy trapped in waveguide Di also decays because part of the energy couples to 
the middle waveguide from which the energy is lost into the environment at a rate y. To avoid 
a strong coupling of energy between Di and D3, we gradually change the distance between the 
two side waveguides and the middle one to guarantee an adiabatic process. In addition, large 
phase mismatchings A13 and A23 are used to suppress the energy coupling to waveguide D3. 
In the output side, we decouple the waveguides Di and D3 by introducing a large distance to 
keep the light energy in Di almost constant. The profile of the mode is also kept stable. The 
couplings used in coupled mode equation for fitting the following numerical results are shown 
in Fig.|2](b). These coupling are strongly dependent on the distance c/i3,23- 

3. Results 

Now we study the left-right nonreciprocity where we have nonreciprocal light transfer in the 
transverse direction ifTS] ISTI . Similar to coherent population trapping in quantum optics, we 
can trap most light energy in the selected waveguide Di by designing a weak coupling in 
comparison with K23. Our numerical results shown in Figs. |3](a) and (c) demonstrate a left- 
right nonreciprocal transverse energy flow. Whatever port Pu or Pl2 we choose to lauch the 
light into, most of the light is trapped in the waveguide Di, and comes out from the same port 
Pri (blue lines). The transmission for light input into port P12 is about 27% (Figs.|3](b)) but it 
increases to 40% if the light is incident into port Pu (Figs.|3](d)). The light in the waveguide D2 
leaks to D3 and subsequently is absorbed as it propagates. The contrast ratios of light intensities 
in waveguides Di and D2 are higher than 29 dB in both cases. 

So far we presented results based on the numerical solution of the Heknholtz equation. Next 
we compare these results with the solution of the coupled mode equation Eq. dUi, see dashed 
lines in Figs. [3] (b) and (d). Clearly, the coupled mode theory agrees well with the numeri- 




Fig. 3. Left-right nonreciprocity corresponding to Fig.|2](b). The field propagates from left 
to right in photonic circuits, (a) Light incident into the waveguide D2; (c) light enters the 
waveguide Di; (b) and (d) Intensities of field at the middle of waveguides Di [blue lines 
(i)], D2 [red lines (ii)] and D3 [green lines (iii)]. Dashed thin lines are the corresponding 
plots by solving Eq. 



cal simulation. However, the coupled mode equations predict power oscillations in the coupled 
regions. Moreover, the intensity of light in waveguide Di in Figs. [3] (b) are higher than the 
numerical results of BPM method. The origin of discrepancy are complex. In our numerical 
simulation, the waveguide is 2 |xm or 4 |j.m wide, and enough to support higher modes. The 
width of waveguides D12 also change sUghtly in the bending region. In this region, the funda- 
mental mode couples to higher modes and lose to environment. There are also many controlling 
parameters in the coupled mode equation Eq. (|4]i. We do not know the exact varying values of 
these parameters corresponding to the structure of system. In addition, we use a simple model, 
Eq. (lU, in which the parameters are different from the values given by a full coupled mode 
theory 1291 f30l . Thus it is difficult to perfectly fit the numerical results with the coupled mode 
theory. In spite of small discrepancy, the coupled mode theory still fits numerical results in 
principle. 

The bandwidth in which the propagating light show nonreciprocal behavior is an important 
feature. We study the frequency dependence of the transmission spectra in Fig.|4]by numerical 
simulations. As the widths and refractive indices of the two side waveguides Di and D2 are the 
same, the propagating constants ySi and ^62 are equal. Thus the phase mismatch An is euqal to 
A23 ideally. As a result, the transmissions are insensitive to the frequency of incident light. Thus 
our system has an ultrabroadband nonreciprocal window over 200 nm, from 1 .6 [a.m to 1 .8 [a.m. 
According to Eqs. ([U and the structure of system is scalable in size to shift this frequency 
window. In the nonreciprocal windows, the light in waveguide D2 is always vanishing because 
it couples to the lossy channel D3. Next we concentrate our discussion in the nonreciprocal 
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Fig. 4. Nonreciprocal transmission as a function of frequency of input light using structure 
as in Fig. [3] (a). Tliin blue lines (i) and (iii) show the light trapped in waveguide Di, thick 
red lines (ii) and (iv) show the light energy in D2. Solid lines for light launching to D2, 
dashed lines for light input into Di . 

window of interest. It can be clearly seen in Fig.|4] whatever waveguide the light is incident to, 
about 30% energy is trapped in D 1 and comes out of port Pri . The weak frequency-dependence 
of transmission comes from the change of eigenmode profiles and their couplings /en and K23, 
which are also dependent on the profiles of eigenmodes and wavelength [ref. to Eq. (|5]l]. The 
deviation in fabrication may result in a small difference between A13 and A23. However the 
transmission change slightly if A 13 » A23. 

4. Conclusion 

In conclusion, using a technique analogous to the coherent population trapping in quantum 
optics, we broke the symmetry of transverse light propagation in the photonic circuits of three 
coupled waveguides. Our proposed system is made only from linear, passive materials. Our 
simulations indicate the possibility of asymmetric transverse energy flow in an ultrabroadband 
window spanning over 200 nm in frequency. Although our proposed system has a relatively 
large insertion loss, it opens a door to the possibility of highly efficient optical nonreciprocity 
in a linear, passive medium. 
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